Density fluctuations on mm and Mpc scales 



N. P. Basse * 

Plasma Science and Fusion Center, Massachusetts Institute of Technology, 

Cambridge, MA 02139, USA 



Abstract 

We will in this paper report on suggestive similarities between density fluctuation 
power versus wavenumber on small (mm) and large (Mpc) scales. 

The small scale measurements were made in fusion plasmas and compared to 
predictions from classical fluid turbulence theory. The data is consistent with the 
dissipative range of 2D turbulence. Alternatively, the results can be fitted to a 
functional form that can not be explained by turbulence theory. 

The large scale measurements were part of the Sloan Digital Sky Survey galaxy 
redshift examination. We found that the equations describing fusion plasmas also 
hold for the galaxy data. 

The comparable dependency of density fluctuation power on wavenumber in fu- 
sion plasmas and galaxies might indicate a common origin of these fluctuations. 
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1 Introduction 



If one were to make a survey of where we are, what we know and what we 
do not know about magnetically confined fusion plasmas, turbulence would 
certainly be an area marked 'Here Be Monsters'. Cross-field transport (per- 
pendicular to the main magnetic field) assuming that only binary particle 
collisions contribute is called neoclassical transport [1]. This transport level 
includes effects associated with toroidal geometry. However, in general the 
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measured transport is several orders of magnitude larger than the neoclassical 
one, especially for electrons. This phenomenon has been dubbed anomalous 
transport and is subject to intense studies on most experimental fusion devices 
[2] . Anomalous transport is beheved to be driven by turbulence in the plasma. 

It is generally thought that turbulence creates fluctuations visible in most 
plasma parameters. Therefore a concerted effort has been devoted to the study 
of fluctuations and their relation to the global (and local) plasma confinement 
quality. 

In this paper we study electron density fluctuation power versus wavenumber 
(also known as the wavenumber spectrum) in the Wendelstein 7-AS (W7- 
AS) stellarator [3]. Wavenumber spectra characterize the nonlinear interaction 
between turbulent modes. The fluctuations were measured using the localized 
turbulence scattering (LOTUS) diagnostic [4,5]. 

As we shall see, the density fluctuation power P decreases exponentially with 
increasing wavenumber k on mm scales in fusion plasmas 

P{k) ocjX e-'^^ (1) 



where n > is a constant having a dimension of length and k = 27r/A, where A 
is the corresponding wavelength. This was initially noted using the simplified 
form 

P{k) oc 6-"^= (2) 



in Ref . [6] . Eq. (2) also holds for density fluctuations in the Tore Supra tokamak 
[7]. 

Having the exponential structure of the wavenumber spectrum in mind, we 
were intrigued to see a flgure in KVANT [8], a magazine published by the 
Danish Physical Society, showing the root-mean-square density fluctuation 
amphtude of galaxies from the Sloan Digital Sky Survey (SDSS) [9] versus 
distance that seemed to display the same exponential behavior as we found 
in fusion plasmas. This remarkable similarity prompted us to apply the same 
analysis to the SDSS data as we had previously used for our fusion plasma 
measurements. 

The paper is organized as follows: In Section 2 we summarize our wavenumber 
spectrum measurements in fusion plasmas. Thereafter we describe inflation 
and the SDSS wavenumber spectrum in Section 3. We discuss the results in 
Section 4 and flnally state our conclusions in Section 5. 
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2 Fusion plasmas 



A wavenumber spectrum of turbulence in W7-AS is shown in Fig. 1. The 
measured points are shown along with two power-law fits 

P{k) oc k-"", (3) 

where m is a dimensionless constant. The power-law fits are shown as solid 
lines and an exponential fit to Eq. (1) is shown as a dashed line. The power- 
law fits are motivated by classical fiuid turbulence theory where one expects 
wavenumber spectra to exhibit power-law behavior with exponents m depend- 
ing on the dimension of the observed turbulence: 

• 3D: Energy is injected at a large scale and redistributed (cascaded) by non- 
linear interactions down to a small dissipative scale. In this case, the energy 
spectrum in the inertial range E{k) oc k~^^^ [10]. 

• 2D: Here, two power-laws exist on either side of the energy injection scale. 
For smaller wavenumbers, the inverse energy cascade obeys E{k) oc k~^^^ 
and for larger wavenumbers, the enstrophy cascade follows E{k) oc k~^ [11]. 

• ID: Energy is injected at a large scale and dissipated at a small scale; 
E{k) oc k-^ [12]. 

Our measured power is equivalent to the ci-dimensional energy spectrum Fd{k) 
[13,10,14] 



P{k) = Fd{k) = 




Ai^2 2TTk A3 - 47rA;^ (4) 

where A^, is the surface area of a sphere having radius k and dimension d. 
Usually on would assume that = 2 in fusion plasmas, since transport along 
magnetic field lines is nearly instantaneous. The fits to Eq. (3) in Fig. 1 yield 
exponents m — 3 (small wavenumbers) and 7 (large wavenumbers) . A similar 
behavior has previously been reported in Ref. [15] where it was speculated that 
the wavenumber value at the transition between the two power-laws should 
correspond to a characteristic spatial scale in the plasma. The only length 
scale close to the transitional value was found to be the ion Larmor radius pi. 



The spectrum at small wavenumbers is roughly consistent with the inverse en- 
ergy cascade in 2D turbulence, F2{k) oc k~^^^. The exponent at large wavenum- 
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Fig. 1. Wavenumber spectrum of broadband turbulence in W7-AS measured using 
the LOTUS diagnostic. Squares are measured points. Solid lines are power-law fits 
to the three smallest and five largest wavenumbers, the dashed line is a fit to Eq. (1). 
The power-law fit grouping of points is the only one where convergence is obtained. 

bers does not fit into this framework. However, for very large wavenumbers 
one enters the dissipation range; here, it has been argued that the energy 
spectrum could have one of the following dependencies 



where a > is a constant having a dimension of length (see Ref. [12] and 
references therein). The energy spectrum proposed by J. von Neumann was 
what initially inspired us to investigate an exponential decay of P{k) in Ref. 
[6]. Fitting all wavenumbers to Eq. (1), -ENeumann(fc)M2, we find that n = 0.1 
cm or a wavenumber of 55 cm~^. Alternatively, the transitional wavenumber 
found at the separation between the two power-laws is 31 cm~^ (0.20 cm). The 
expression £^Heisenberg(fc)M2 yiclds m = 8, which is close to the experimental 
value m = 7 for large wavenumbers. Calculating the ion Larmor radius at the 
electron temperature ps for this case we find that it is 0.1 cm, i.e. the same 
order of magnitude as the spatial scales found above. We used ps instead of pi 
because ion temperature measurements were unavailable. 

Currently we can think of three possible explanations for the behavior of the 
wavenumber spectrum: 

(1) We observe 2D turbulence and the transition between the two power-laws 
occurs at a spatial scale where the inverse energy cascade develops into 
the dissipation range. However, the enstrophy cascade is not accounted 



-^Neumann 

(k) oc e 



—ak 




(5) 
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for in this case. 

(2) We observe 2D turbulence in the dissipation range described by a single 
exponential function as proposed by J. von Neumann. 

(3) Turbulence theory does not apply. The transition between two power- 
laws or the characteristic scale found using a single exponential function 
(Eq. (2)) indicates that one scale dominates the turbulent dynamics in 
the wavenumber range studied. 



3 GalELxies 

3. 1 Inflation 

Dramatic developments have taken place in cosmology over the last decade, 
lending increasing support to the paradigm of inflation as an explanation 
for what took place before the events described by the big bang theory [16]. 
Inflation solved the so-called horizon and flatness problems, but was at odds 
with earlier observations indicating that the ratio of the mass density of the 
universe to the critical value, the density parameter Q, was 0.2-0.3, while 
inflation predicted it should be 1: 

Pc 

Q, <1 : open 0,-1: flat Q, > 1 : closed, (6) 

where pc — SHq/SttG is the critical mass density, Hq — 70 km/s/Mpc is the 
Hubble parameter observed today and G is I. Newton's gravitational con- 
stant [17]. However, new measurements in the late 1990's lead to a drastic 
modification of Q: Observations of type la supernovae (SN) showed that the 
separation velocity between galaxies was speeding up, not slowing down as 
would be expected for an open universe. The underlying explanation for this 
accelerated expansion is not understood, but it seems that the universe con- 
tains large quantities of negative pressure substance, creating a gravitational 
repulsion driving the expansion. This negative pressure material is called dark 
energy, the total density of dark energy fl\ is 0.7. The existence of dark en- 
ergy is equivalent to the cosmological constant A introduced by A. Einstein. 
The dark matter density is 0.25 and the baryonic matter density D,h is 
0.05, so the total density is very close (or equal) to the critical density. The 
SN la data is supported by measurements of nonuniformities in the cosmic 
microwave background (CMB) radiation. The CMB anisotropy is due to the 
presence of tiny primordial density fluctuations at the time of recombination. 
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Fig. 2. Wavenumber spectrum of galaxies measured by the SDSS Team. Squares 
are measured points. Solid lines are power-law fits, the dashed line is a fit to Eq. 
(1). The power-law fit grouping of points is chosen so the combined, normalized, 
of the fits is minimized. The data is taken from Ref. [23]. 

where atoms formed. At that point in time the age of the universe was about 
300,000 years and the temperature was 3000 K. The structures observed in 
the CMB are called acoustic peaks, and the simplest versions of inflation all 
reproduce these structures quite accurately. The acoustic peaks can not be 
modelled by assuming that the universe is open. 



3.2 Wavenumber spectrum 

A study of density fluctuations on large scales using 205,443 galaxies has been 
published by the SDSS Team in Ref. [18], see Fig. 2. 3D maps of the universe 
are provided by the SDSS galaxy redshift survey, observing about a quarter of 
the celestial sphere using a 2.5 m telescope and a charge-coupled device (CCD) 
camera. The galaxies had a mean redshift z ~ 0.1, corresponding to light 
emitted 1-2 Gyr ago [17]. Fixing some cosmological parameters to Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite values [19,20,21] one finds - 
using physics based models - that the wavenumber spectrum measurements 
were fitted by a matter density flm = + = 0.295 ± 0.0323. In this case 
h = Hq/[100 km/s/Mpc) = 0.72 was assumed and it was observed that the 
wavenumber spectrum was not well characterized by a single power-law. 
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A follow-up paper by the SDSS Team, Ref. [22], combined non-CMB mea- 
surements (SDSS) with CMB measurements (WMAP) to constrain free pa- 
rameters of cosmological models and break CMB degeneracies in parameter 
space. This resulted in = 0.30 ± 0.04 and h = 0.70l|]:g;^. Adding the SDSS 
information more than halved WMAP-only error bars on some parameters, 
e.g. the Hubble parameter and matter density. 

The data presented in Fig. 2 has been taken from M. Tegmark's homepage 
[23]. According to the recommendation by the SDSS Team [18], the three 
largest wavenumbers shown are not used in the fits described below. 

As we did for the W7-AS data in Section 2, we fit the SDSS measurements 
to two power-laws (Eq. (3)) or a single exponential function (Eq. (1)). The 
power-law fits are shown as solid lines, the exponential fit as a dashed line. 

The power-law fits yield exponents m = 0.8 (small wavenumbers) and 1.4 
(large wavenumbers). The wavenumber ranges were determined by minimizing 
the combined, normalized, of the fits. As the SDSS Team found, a single 
power-law can not describe the observations. The exponents are not close to 
the ones governing fiuid turbulence discussed in Section 2. The transitional 
wavenumber is 0.09 h Mpc~^, corresponding to a length of 67 h~^ Mpc. 

We find the characteristic length from an exponential fit to be n = 2 h~^ Mpc 
or a wavenumber of 3 h Mpc~^. 



4 Discussion 

The fact that density fiuctuations on small (fusion plasmas) and large (galax- 
ies) scales can be described by an exponential function might indicate that 
plasma turbulence at early times has been expanded to cosmological propor- 
tions. A natural consequence of that thought would be to investigate fluctu- 
ations in quark-gluon plasmas (QGPs) corresponding to even earlier times. 
However, experimental techniques to do this are not sufficiently developed at 
the moment due to the extreme nature of QGPs. 

It is fascinating that wavenumber spectra over wider scales peak at small 
wavenumbers and decrease both above and below that peak. This is seen both 
in fusion plasmas [24] and for galaxies, see e.g. Fig. 38 in Ref. [18]. Turbulence 
theory in ID or 3D would interpret the peak position as the scale where energy 
is injected. 

Fitting wavenumber spectra to power-laws is based on fluid turbulence the- 
ories, but in general care must be taken when interpreting the outcome: We 
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know that an exponential function can be Taylor expanded to an infinite power 
series: 



So locally, i.e. for a small range of wavenumbers, an exponential dependency 
can be masked as a power-law; the exponent would vary as a function of the 
wavenumber range selected. 

We favor the exponential functions over power-laws as descriptors of the data, 
either including an algebraic prefactor as in Eq. (1) or using a pure exponential 
function as in Eq. (2): In the latter case, fits yield n = 0.2 cm for fusion plasmas 
and 11 h~^ Mpc for galaxies. The exponential decrease of density fluctuation 
power versus wavenumber implies that either (2) or (3) in Section 2 could 
explain the data. Perhaps forcing occurs at a large scale and transitions either 
directly to dissipation (2) or to some effect governed by the scale determined 
from the exponential function (3) at smaller scales. The fact that we obtain 
n ^ Ps using Eq. (1) for fusion plasmas supports this conjecture. 



5 Conclusions 

We have in this paper reported on suggestive similarities between density 
fluctuation power versus wavenumber on small (mm) and large (Mpc) scales. 

The small scale measurements were made in fusion plasmas and compared to 
predictions from turbulence theory. The data flt Eq. (1), which is consistent 
with the dissipative range of 2D turbulence. Alternatively, the results flt Eq. 
(2) which has a functional form that can not be explained by turbulence theory. 

The large scale measurements were part of the SDSS galaxy redshift survey 
As is the case for fusion plasmas, the galaxy data can be described by Eqs. (1) 
and (2). The similar dependency of density fluctuation power on wavenumber 
might indicate a common origin of these fluctuations, perhaps from fluctua- 
tions in QGPs at early stages in the formation of the universe. 

The cross-disciplinary work presented here is hopefully just the beginning of 
an interesting path that can bcneflt both flelds. As a flrst step, we will expand 
our studies to encompass a wider range of scales, both for fusion plasma and 
galaxy measurements. 




i=0 



(7) 
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